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The microscopic origin of triplet superconductivity in Sr 2 RuC>4 is discussed, paying atten- 
tion to the role of Coulomb interaction, U pp , at the O site. It is shown on the d-p model that 
U pp induces a ferromagnetic exchange interaction between " d-electrons" (molecular orbital 
with dry-symmetry) at adjacent Ru sites, leading to short-range ferromagnetic correlations 
and promoting Cooper pairing with (sinf> x ± i sinpj / )-symmetry on the 7-band. The reason 
why such ferromagnetic correlations work effectively may be traced back to the fact that the 
level of 4<i-electrons at Ru sites is relatively low and located near that of 2p-electrons at O 
sites. 
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Since its discovery by Maeno et al, superconductivity in Sr2Ru04 has attracted much 
attention both experimentally and theoretically. Now it seems to have been established ex- 
perimentally that the superconducting (SC) state is in the chiral spin triplet state. 2,3 The 
existence of technically linelike nodes in the SC gap has also been suggested. 4-6 Since these 
observations are in contrast with the expectation at the early stage of theoretical works, ' 8 
some phenomenological gap models have been proposed to explain those experiments. 9-11 The 
results of ref. 9 were consistent with the observed gap with linelike nodes in [100] and [010] 
directions. 12 On the other hand, microscopic calculations have been performed on various 
models and methods. 13-16 The results on the multiband Hubbard model by Nomura and Ya- 
mada 15 seem consistent with almost all the available measurements to date probing the gap 
structure. However, the Hubbard model with only on-site repulsion, U^, seems too stoical 
to discuss the superconductivity in Sr2Ru04 because the quasi-particles (QP) consist of 'id- 
electrons at Ru sites and 2p-ones at O sites. Then the intersite interaction between QP, arising 
through the Coulomb interaction at O sites, may not be neglected, while the direct Coulomb 
interaction V between 4cf-electrons at the nearest-neighbor Ru sites 17 would be negligibly 
small. The interaction overlooked so far is 
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where p m %a is the operator for a 2p m -electron (m = x, y) of spin a at the ith site, and <7 Q( 3's are 
the Pauli matrices. Such a Coulomb interaction gives rise to an effective intersite correlation 
J between the molecular orbital with d^-symmetry ("d-electron"). We expect that eq. (1) 
can be expressed phenomenologically as 

W cx ^ -^2 ^J^+^^^d^sd^ia^-a^s), (2) 

kk'q a/3-yS 

where d denotes the operator for a "d-electron" and J q = 2J(cosq x + cosq y ). Although 
exchange interaction (2) is a too simplified version, we derive a more realistic and complicated 
version later starting with the d-p model and show that this kind of intersite interaction 
promotes short-range ferromagnetic correlations (SRFMCs), which have been measured quite 
recently by inelastic neutron scattering. 18 The purpose of this Letter is to demonstrate that 
SRFMC promotes triplet pairing with (sinp x ± isinp 3/ )-symmetry, as proposed in ref. 9. We 
follow the weak-coupling calculations with the pairing interaction given by the second-order 
perturbation theory (SOPT) with respect to Udd and U pp . 

Let us consider the situation drawn in Fig. 1, where a 2p x - or 2p^-orbital is sandwiched 
between two d-orbitals with d^-symmetry. Hereafter, we restrict our discussion within the 




Fig. 1. Schematic figure of Ru02 plane. td P and t pp are hopping integrals. Gray domains circled by 
dashed squares represent the molecular orbitals, which form a two-dimensional square lattice. 



7-band since it is expected to provide the main contribution. 15 The noninteracting d-p Hamil- 
tonian, which can well describe the situation shown in Fig. 1, is given as follows: 
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where d\^ a and p m \to are the operators for 4d- and 2p m -electrons of momentum k and spin 
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a, respectively, V m y_ = 2ifd p sin(^p-), and = 4t pp sin(4f) sin(^-). Note that we should 
take t pp < for discussing Sr2Ru04, while t pp > for the cuprates. Solving the equation 
det^ - A k i) = 0, we can obtain the eigenvalue Ajk and eigenvector (i = 1, 2, 3). Then, 
the diagonalized Hamiltonian becomes 

H? = EW, (4) 

k,cr 

Here, </>ko- is related to V'ko- by the relation (f)^ = U^ 1 ^^, where the matrix is defined as 
U\s_ = («ik «2k ^3k)- Substituting the relations between 0k CT and ipk a into eq. (1), we obtain, 
as shown below, the effective interaction between QP described by a^. For this purpose, it 
is enough to confine our discussions within the case t pp = 0: namely we need not consider 
the hopping t pp . Because electrons located on both sides of a certain O site can interact 
with each other effectively only through the O site, the term including t pp corresponds to 
the next nearest-neighbor interaction (higher-order contribution). Then, hereafter, we can use 
simplified relations such as 

, Aik A 2 k . A 3k , , 

«kcr = 7 Olkcr + 7 02ko- + 7 03ko-, (OJ 

-^lk -^2k -^3k 

^Tnk ^Trak ^Tnk / a \ 

Pmka ~ 7 a lka + "T a 2ka + 7 a 3ka, (P) 

-^lk -^2k ^3k 

where = (iV^k) 2 + |K/k| 2 + -^fk) 1 ^ 2 (* = 1> ^> Moreover, it is sufficient to consider only 

the QP on the orbital i = 1, because the other bands are completely filled with electrons. 

Hereafter, we denote aik CT and Lik simply as a^ a and Lk, respectively. The operator a^ a is 

regarded as the dko- i n eq. (2). 

Now we introduce the interaction terms. Using relations (5) and (6), we express such terms 

using ako-'s. Substituting relation (6), into the Fourier transformed form of eq. (1), we obtain 

u(x) U pp sr^ ^rk+qKck'-qKck'Krk f + , 

nU = -yLL L k+q L k , q L k ,L k < + «A'-^saw(v a p-a, 5 ). (7) 

At first sight, it seems that this interaction Hamiltonian, eq. (7), has a complicated form in 
real space. However, if we neglect the k-dependence of Lk's where Lk = A ~ 0(ed — e p ), we 
can perform the Fourier transformation of eq. (7), leading to 

Hcx = 3A 4 ^ ^ S 

+ (i,3,i,j) + (hhjj) ~ ~ (hhj,i) 

- ti,i,i,i) ~ (i,3,i,i)}(vai3-v~ts), (8) 

where we have added the term arising from the exchange process along the y-direction, the sum 
is taken within the nearest-neighbor sites, and the abbreviation k, I) represents the 
term a^a^aksai^. Eventually, it is revealed that eq. (7) contains all the types of interaction 
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between two electrons at nearest-neighbor Ru sites. Performing the Fourier transformation 
again and executing straightforward calculations, we can rewrite eq. (8) into the following 
form: 

W ex = — - ^2 ^k,k';q a k +qQ «k'_ q7 ak'<5«k/3(o ; a / 3-0 ; 75 ~ ^a^-yd), (9) 

kk'q a/?7<5 

J k) k' ;q = 7 z PP dp T — p- Y]\ - + cosq m - cos k m - cos k' m + - cos(A; m + k m )\. (10) 

Mc+q^k'-q-kk'^k ^ L ^ ^ J 

In eq. (10), the k-dependence of £k's is recovered again. The second term in the curly brackets 
of eq. (10), cos ({mi corresponds to the simplified version of eq. (2). There are many terms which 
are not included in the naive estimation, eq. (2). Therefore, the coupling constant J cannot 
be written simply as J q , but has a complicated momentum dependence Jk,k';q as eq. (10). 
The term S a p5-ys m eq. (9), which appears automatically through the procedure above, assures 
that, due to the Pauli principle, there is no interaction between electrons with the parallel 
spin direction originally on the same O site. For interaction (9) to play an important role, the 




ratio tdp/(£d ~ £p) = tdpl 'A should not be too small (namely, A/td p should not be too large). 
Such a condition is expected to be fulfilled in Ru-oxides because the level of 4d-electrons 
in Ru is deeper than that of 3d-electrons in Fe or of 5(f-electrons in Os in general, and the 
level is located near that of 2j>-electrons at the O site. This implies that there exists a strong 
hybridization between Ru and O. Indeed, as shown in Fig. 2, the partial density of states 
(DOS) for 4:d xy , 2p x and 2p y obtained by band structure calculations 19,20 are reproduced by 
taking A/td p ^ 0.15 with the use of eqs. (5) and (6): namely, p 7 ~ 2.0po (Py — 0.5pR u ) where 
Py> PRu and po denote the partial DOS per atom (at the Fermi level) of Ru (7-band), Ru 
(total of a-, (3-, and 7-bands) and O (in-plane), respectively. The same procedure is applicable 
to the on-site Coulomb repulsion Udd between ci-electrons, Udd^id\^d^d\,di^. Along the way 
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parallel to obtaining exchange interaction (7) from the original interaction between p-electrons, 
eq. (1), the Hubbard interaction is expressed in terms of the operators describing the QP 
band. Then, the total Hamiltonian is given by 

ft = X]( Ak ~ ^)«L a k<x + ^2 ^k,k';q4+qT a k'-q| a k'l a kT + ^ex, (H) 
k.cr kk'q 



^k,k';q = — j r t ' \ 12j > 



^dAk+qAk'-qAk'A k 
8-^k+q^k' -q^k' -^k 

where [i denotes the chemical potential. It is noted that the repulsion U^d m eq- (12) should 
be regarded as an effective interaction between QP, that is reduced by the electron correlation 
leaving a relatively heavy QP. We set t pp = —OAtdp to reproduce the Fermi surface of the 
7-band given by the band structure calculations. Then the bandwidth of the 7-band is W ~ 
2.2tdp. The last two terms of eq. (11) can be rewritten into a more convenient form to perform 
the perturbation expansion as follows: 

Hint = ^k,k'; q a k +qT a k'-q| a k'l a kT, (13) 
kk'q 

^k,k';q = £4:,k';q + ^k,k';q + ^k',k;k-k'+q- (14) 

Hereafter, td p is adopted unit of energy. 

The static spin susceptibility x_i_(k, 0) including vertex correction, up to the first-order 
perturbation, is given by 

X±(q, kj) = -T^jl-T^ Jp + q,p';-p + p'S(p' + q, ie / +iw)0(p',te)}0(p + q, ie+iw)g(p, is), 
£,P e',p' 

(15) 

whose diagrammatic expression is shown in Fig. 3(b). It is enough to calculate x_L because 
the SU(2) symmetry is preserved. The results of numerical calculations are shown in Fig. 4. 



(a) (b) 

Fig. 3. (a) Vertex Jk,k':q- (b) Spin susceptibility X_i_(q, iu>). 



X±(q, 0) has some peaks at some incommensurate wave vectors Qi c 's. The Fermi surface 
simulating the 7-band of Si'2Ru04 has no nesting tendency, although the Fermi surface shifted 
by Qi c 's has a point contact in the extended Brillouin zone scheme. 21 Therefore, the peak 
structure of XJ-(q>0) at Qic's is considered to come from the band effect. In particular, the 
shape of X_!-(q>0) consists of a broad peak around q = (0,0) and a small dip (at q ~ (0,0)), 
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which gives rise to a practical SRFMC. The effect of Udd is to enhance X_i_(q, 0) in the entire q- 
space as seen in Fig. 4. On the other hand, the peak around q = (0, 0) is much more enhanced 
by U pp than by Udd, while there is no extra increase by U pp at (ir, tt). The qualitative behavior 
of x±(Q)0) remains unchanged if we neglect the k-dependence of L^'s in eq. (10) or (12). 
Then, we can draw the picture that SRFMC is induced by the effect of on-site repulsion 
between p-electrons at the O site. The broad peak structure at q = (0, 0) is consistent with 




Fig. 4. Momentum dependence of static spin susceptibility XJ-(q, 0) for several values of Udd and U pp . 
The dashed lines that have no parameter values correspond to, from the lower line, (Udd, U pp ) = 
(1, 0), (2, 0) and so on, respectively. Namely, we first applied Udd, then applied U pp . 



the SRFMC that has been measured quite recently by inelastic neutron scattering, 18 and is 
responsible for the occurrence of triplet superconductivity in Sr2Ru04 as discussed below. As 
for the results of different electron fillings, the more electrons are doped, the more the peak 
structure around q = (0, 0) is enhanced. 

Diagrams for the irreducible pairing vertices up to SOPT in terms of J are shown in Fig. 5. 
Its analytic expressions are given as follows: 



•t t- 



t t 



= \ ■ 
•t t-» 



V 



r 




Fig. 5. Pairing interaction vertices and up to the second-order perturbation with respect to 
J. The solid line and the square denote the noninteracting Matsubara Green function Q and the 
intcrsitc coupling vertex J, respectively. 



r k|k' = ^k,p;-k+k'«/ P +k-k' ,-k ; -k+k'£(p, ie)£(p + k - k', ie), (16) 

e,p 
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r k|k' = ^k ,-k;-k+k' - Jk,p;p+k'</-k,p+k+k';p+k£(P, is)G(p + k + k', fc). (17) 

e,p 

The pairing interaction in the triplet manifold, V^ k ,, is given by eq. (16) directly, and that in 
the singlet manifold, V^ k ,, is given by 

^k,k' = Jk -k;-k+k' - <4,p;p-k'-/-k,p+k-k';p+k£(P, '^)G(p + k - k', ifi). (18) 

e,p 

Note that SU(2) symmetry is preserved in the triplet manifold. Namely, the triplet component 
stemming from eq. (17) is the same as the one from eq. (16). 

At T = T c , we can use the linearized gap equation, in the weak-coupling formalism, given 

by 

A k = -£^t^(i)' (19) 

A proper way of obtaining the gap structure with the highest T c is to solve eq. (19) as an 
eigenvalue problem without specifying the type of the gap A^. However, for simplicity, we 
seek the variational solution for eq. (19). Namely, we calculate T c by specifying the type of 
gap function, such as \/2sin/c x (p^-pairing) , y/2 sm(k x + k y ) (pz+y-pairing) for triplet pairing, 
and 2 sin k x sin k y (d^-pairing) , [cos k x — cosfe,,] [d x i_ y i -pairing) for singlet pairing. However, 
it is rational to consider that the type of pairing giving the highest T c is dominant in the 
true gap function. The pairing interaction V^ k ,, eq. (16), and V^,, eq. (18) is estimated at 
T = 0.02, because their values are not sensitive at temperatures T < 0.02. We adopt the first 
Brillouin zone divided into 45x45 k- meshes. The results for a series of parameters are shown 
in Fig. 6. We can see in Fig. 6 that T c is enhanced as U pp is applied. For a system with a 
fillings of n = 1.33, a small value of A (splitting between p- and d-levels), and a moderate 
value of Upp, corresponding to Sr2Ru04, p^-pairing would be realized as shown in Fig. 6(a). 
When the system is located away from the half-filling state, it is also enhanced as shown in 
Fig. 6(b). On the other hand, p x+y -paivmg has T c less than 10~ 3 . In ref. 17, it is discussed 
that the next nearest-neighbor (diagonal) pairings, d xy - or p^+^-pairing, are promoted due to 
the nearest-neighbor Coulomb repulsion V . In our case, on the contrary, the nearest-neighbor 
pairing, p x -pairing, is realized through SRFMC of H cx . As for the other symmetries, d x 2_ y 2 
and d xy , we cannot obtain a finite T c in the parameter region shown in Fig. 6. At first sight, 
it appears unrealistic that U pp is larger than Udd ~ W, because 3<i-electrons are much more 
localized on the ion site than 2p-electrons. However, it is crucial to note that the Coulomb 
interactions here are effective ones that are greatly renormalized. Indeed, the effect of U pp on 
QP cannot be screened by avoiding U pp because QP consisting of molecular orbitals at the 
nearest-neighbor sites cannot get rid of U pp . On the other hand, the effective interaction Udd 
can be reduced by correlations and by making a heavy QP. Therefore the values adopted for 
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Fig. 6. T c for triplet pairing with p^-symmetry as a function of U pp . (a) A dependence, (b) doping 
dependence. 



Udd and Upp are considered to be not unrealistic. 

In conclusion, we have derived from the so-called d-p model the effective intersite exchange 
interaction TC ex , whose origin is the on-site Coulomb interaction at an O site. SRFMC were 
induced by this exchange interaction. Then, we have shown that the triplet SC state of (s'mp x ± 
isinp 3/ )-symmetry is promoted as applying U pp within SOPT. These results simulate quite well 
the properties observed in Sr2Ru04. While our results are derived on SOPT for the pairing 
interaction, they might make sense at least in describing the qualitative behavior. 22 It is 
possible to take into account the effect of a- and /3-bands on the same footing as the 7-band 
since the exchange interaction arising from the same mechanism also works between QPs in 
these bands. According to preliminary calculations, the pairing with [sin(p x +p y ) +isin(p x — 
p^)]-symmetry could be promoted, which is consistent with a line-node-like gap along the 
diagonal direction k x = ±k y . 12,15 
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